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We report electrical tuning by the Stark effect of the excited-state structure of single nitrogen- 
vacancy (NV) centers located < 100 nm from the diamond surface. The zero-phonon line (ZPL) 
emission frequency is controUably varied over a range of 300 GHz. Using high-resolution emission 
spectroscopy, we observe electrical tuning of the strengths of both cycling and spin-altering transi- 
tions. Under resonant excitation, we apply dynamic feedback to stabilize the ZPL frequency. The 
transition is locked over several minutes and drifts of the peak position on timescales ^100 ms are 
reduced to a fraction of the single-scan linewidth, with standard deviation as low as 16 MHz (ob- 
tained for an NV in bulk, ultra-pure diamond). These techniques should improve the entanglement 
success probability in quantum communications protocols. 



Integrated photonic networks based on cavity- 
coupled solid-state spin impurities offer a promising 
platform for scalable quantum computing [T]-[5]. A 
key ingredient for this technology is the generation 
and interference of indistinguishable photons emitted 
by pairs of identical spin qubits [BHS]- This requires 
spectrally stable emitters with identical level structure, 
a formidable challenge in the solid-state environment. 

A potential solution is to use external control to coun- 
teract sample inhomogeneities. In candidate systems 
based on single molecules [9HTT], quantum dots [12^ 13j, 
and negatively-charged nitrogen-vacancy (NV) centers 
in diamond [ 14HT6] . the level structure can be statically 
tuned via the DC Stark effect. However, the spectral 
stability of emitters in these systems is often hampered 
by local fluctuations which cause the emission frequency 
to change with time, a phenomenon known as spec- 
tral diffusion [17J. Previous attempts to address this 
problem have focused on improving the host material 
[m [T8H20J or using post-selection techniques [16, |2TI - 
[23] . but a robust, high-yield solution is still lacking. 

The diamond NV center is an attractive spin qubit, 
as it exhibits a unique combination of long-lived spin 
coherence pij and efficient optical control and readout 
[25| [26j . However integration into on-chip photonic net- 
works requires NV centers to be located near nanostruc- 
tured surfaces, where inhomogenous strain and spectral 
diffusion can be particularly problematic |27l |28j. In 
this Letter, we first demonstrate electrical control over 
the zero-phonon line (ZPL) transition frequencies, as 
well as probabilities for both cycling and A-type tran- 
sitions, of single NV centers located near the diamond 
surface. We then show that spectral diffusion of the 
ZPL can be suppressed to 16 MHz standard deviation, 
on timescales ^100 ms, by providing rapid electrical 
feedback to compensate for local field fluctuations. 

The negatively-charged nitrogen- vacancy (NV) center 
has C^^v symmetry, and the basic energy structure is de- 
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FIG. 1. (a) NV level structure at high transverse field 
{d±F± = S± > 15 GHz). Dashed lines indicate the spin- 
selective decay path responsible for optical pumping, (b) 
Effect of longitudinal and (c) transverse electric fields on 
the excited state levels, (d) Fluorescence micrograph of the 
electrode structure. NVl was illuminated by 532 nm light 
and appears white, while the position of NV2 is denoted 
with an asterisk. Metal electrodes appear as shadows. 



picted in Fig. [TJa). The spin-triplet ground state, ^^2, 
is split such that the = spin projection (labeled |0) 
throughout) is separated from the degenerate = ±1 
levels (I ± 1)) by Dqs = 2.88 GHz at T < 100 K [29]. 
The optically excited state (ES) is a spin triplet and or- 
bital doublet with symmetry, and its fine structure 
has been studied theoretically in detail [30, ^J. The 
Hamiltonian describing the ES manifold is: 



JCeS = ^y^SO + ^SS + ^Stark, 
^Stark = —d ■ 



(1) 



where J^so^ ^ssi and J^stark contain respectively, the 
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spin-orbit, spin-spin, and Stark effect contributions; d 
is the electric dipole moment, and F is the electric field. 
The effect of a strain is treated as an effective electric 
field [31 [33]. 

We first consider the influence of J^stark on only the 
orbital portion of the ES wavefunction, consisting of 
two eigenstates, {\(ox)^ Wy})^ initially degenerate at zero 
field. Under electric fields, the orbitals exhibit energy 
shifts, A(F), of: 



A^jF)=d||F||+dxFx, 
A.{F)=duFn-d^F^, 



(2) 



where the directions are with respect to the NV sym- 
metry axis. Longitudinal fields do not lift the orbital 
degeneracy and result only in equal, linear shifts of 
all levels. Transverse fields split the orbitals into two 
branches with an energy difference, 26 ± = (A^^ — A^^), 
that grows linearly with increasing field. The spacings 
between ground-state sublevels remain relatively unaf- 
fected by electric fields [341 135] . The ground state may 
have a longitudinal dipole moment, dcs^W [30) \3T\, but 
in experiments we only resolved Ad\\ = d\\ — dGs,\\^ 

Incorporating spin interactions results in a set of six 
eigenstates, {\A2),\Ai),\Ea:),\Ey),\Ei),\E2)}, ordered 
from highest to lowest energy (at low field). Figures 
[TJb) and (c) show the effect of J^es on all six ES ener- 
gies due to, respectively, Fy and F±. 

We focused most of our study on NV centers close 
to the diamond surface, a necessary feature for future 
integration with nanophotonic devices. Our sample, de- 
scribed in detail elsewhere [36] |33 , consists of a high- 
purity single-crystal, [100]-oriented diamond substrate 
with a ^ 100 nm thick chemical- vapor-deposition-grown 
layer with [NV]?^ 10^ cm"^. The two NV centers stud- 
ied in this work, labeled NVl and NV2, are located in 
this surface layer [37J. Lithographically-defined metal 
electrodes [15j were deposited on the surface [Fig. [TJd)]. 
The layout of the electrodes (labeled Vi,V2 and Vref) 
permits tuning of electric fields in any in-plane direc- 
tion near the center of the structure [Supplementary 
Information] . 

A confocal microscope was used to excite and collect 
emitted light from a diamond sample in thermal con- 
tact with the cold finger of a flow-through, liquid-helium 
cryostat. The cold finger was maintained at a tempera- 
ture T ^ 7 K, and no magnetic field was applied. Two 
forms of spectroscopy were employed: high-resolution 
emission spectroscopy and photoluminescence excita- 
tion spectroscopy. 

For emission spectroscopy, 2 mW of 532-nm laser 
light was focused by a 0.6- numerical- aperture objective 
onto NVl, exciting through the phonon sideband (PSB) 
near saturation. The collected emission was spectrally 
filtered to direct ZPL light (636-638 nm) to a high- 
resolution grating spectrometer. The optical polariza- 
tion was chosen to ensure excitation of both orbital 
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FIG. 2. (a) Stark emission spectroscopy. Va was slowly 
scanned (1 V/min) and emission spectra were obtained in 
one-minute intervals. A dark exposure frame was taken ev- 
ery 8th frame (cyan vertical stripes). Emission frequencies 
are relative to 470.45 THz (637.25 nm). Low-field data are 
inset, (b) Emission spectrum for = 1 V with expected 
peak positions labeled. The \Ex,y) |0) emission line is fit 
with a Gaussian profile, (c) Low-field peak positions from 
(a) and global fit based on Eq. ([T]). Lorentzian fit uncer- 
tainty is smaller than the plotted symbols, (d) Total ZPL 
emission versus 5±. The emission rate was calculated by 
subtracting the mean background and summing counts over 
a range of ~ 3 FWHM linewidths centered at each peak. 
Error bars are based on Poissonian noise, (e) Relative in- 
tensity of the \(oy) ^ emission line along with fit (see 
text). 



branches [38j. Figure [2| a) shows the emission spectra 
versus voltage, Va^ applied simultaneously on Vi and 
Fref, with V2 = 0. By varying Va from to 150 V, we 
observe linear tuning of emission lines over a range ex- 
ceeding 300 GHz. Such a wide tuning range, enabled 
by the enhanced fields provided by our devices [Sup- 
plementary Information], is essential to compensate for 
the large intrinsic fields typical in nanophotonic devices 
[28J. 

Depending on the applied voltage, we resolve between 
one and three emission lines. At K ^ 1 V we observe a 
single emission line [Fig. ^h)] with full-width-at-half- 
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maximum (FWHM) linewidth of 1.4(2) GHz, near the 
spectrometer resolution of ~ 0.9 GHz. We interpret 
this peak as containing unresolved contributions from 
the \Ex^y) |0) cycling transitions [311 HOI- Taking 
into consideration the absence of other peaks, in partic- 
ular the 1^42) I i 1) cycling transition [23j, and the 
observed noise floor, we place a bound on the ground- 
state spin polarization ^^gs = ^0/(^0 + ^±1) ^ 90%, 
where Pi is the occupation probability of state 

Upon application of transverse fields, the spin 
character of the levels in the upper l^x) orbital 
branch, {| A2), | Ai), remain relatively unper- 

turbed. In contrast, the spin character of levels in \S'y)^ 
{\Ey)^ \Ei)^ 1^2)}, niix at avoided crossings due to spin- 
spin interaction [3Ql|3T|[39 , 41j, making these levels use- 
ful for spin-altering A schemes. 

In the range 5 V < 14 ^ 35 V, three emission 
lines are visible [inset of Fig. [2|a)]. Based on the 
positive, linear tuning, the upper peak is identified as 
\Ex) |0). Lorentzian fits to these spectra reveal 
that the two lowest lines are on average separated by 
2.7(2) GHz, which is comparable to Dgs- Consider- 
ing also the negative, linear tuning, we conclude that 
these peaks arise from \S'y) |0) and \(oy) I i 1) 
emission (the three levels within IS'y) are nearly de- 
generate and unresolved here). The presence of these 
lines was previously predicted based on observations of 
spin-altering A-type transitions involving the lower or- 
bital branch [3l SH |42]. Figure j^Jc) plots the emis- 
sion frequencies along with a fit using a model based 
on Eq. ([1]), showing excellent agreement. The fit- 
ted parameters are Ad\\F\\/Va = 0.42(2) GHz/V and 
d±F^/Va = 1.03(3) GHz/V. 

Even with significant emission to | ± 1), we still do 
not observe IA2) I i 1) emission. Throughout, we 
find ^GS ^ 85%. A likely explanation is that any pop- 
ulation in I ± 1) is quickly transferred to the metastable 
singlet levels [43j, preventing the detection of = ±1 
emission lines. This is consistent with Fig. [2|d), where 
the total ZPL emission rate integrated over all lines is 
plotted as a function of one half the orbital splitting. 
Between 3 < S± < 10 GHz (3 < K < 10 V), 
the emission rate falls precipitously before leveling off 
at less than half the initial rate. 

The relative intensity of the emission lines gives fur- 
ther insight into the ES properties. Figure [2|e) plots 
the intensity of the \(oy) | i 1) emission line, normal- 
ized by the total emission from as a function of 5±. 
Evidently, the applied field is a powerful knob in tuning 
the relative transition strengths in this A system. Two 
peaks for the emission of \(oy) | ± 1) are present at 
S_\_ ^ 7 and 15 GHz. These features correspond to level 
anticrossings [see Fig. [ijc)], where maximal mixing of 
levels in the lower orbital branch occurs. The degree 
of mixing depends sensitively on both the magnitude of 
the transverse electric field and its angle, Or^ with re- 
spect to the C^v refiection planes ^39j. We model the 



relative emission intensity by assuming the NV center 
is excited from |0) to one of the three levels in the lower 
branch, \(oy^i). The probability that emission is back 
to I ± 1) is then Po,i{Or, - Po,i{Or, S±)), where 

Po,i{^r^^±) is calculated by taking the overlap of \(^y^i) 
with rus = and tracing over orbital degrees of free- 
dom. Here we assume all levels in \(oy) couple equally 
to the singlets. Using the model based on Eq. ([T]), we 
fit this formula to the data and find good agreement for 
Or = 15(5)°. 

Finally, we also observe a strong dependence of the 
relative emission between the upper and lower branches 
on 5± [Fig. [2|a)]. Given the low temperature T = 7 K, 
this may be due to a single-phonon orbital relaxation 
process. Phonon decay to the lower branch could con- 
tribute to the decreased total emission in Fig. [2|d). A 
detailed study will be the focus of future work. All of 
the effects described above were reproduced in subse- 
quent voltage scans [Supplementary Information]. 

To realize even higher spectral resolution, we 
performed photoluminescence excitation (PLE) spec- 
troscopy. Attenuated light (~ 60 nW) from a tunable, 
external-cavity diode laser (~ 637 nm) was used for ZPL 
excitation near saturation, and the collected light was 
filtered to direct PSB emission (650-800 nm) to a single- 
photon-counting detector. Microwaves resonant with 
the ground-state spin transition, Dgs = 2.877 GHz, 
were continuously applied to counteract optical pump- 
ing [39l [44 j . and light from a repump laser (532 nm) 
was occasionally employed to reverse photoionization 

Figure [sja) plots PLE spectra for NVl as a func- 
tion of Va^ applied simultaneously to V^ei and Vi, with 
¥2=0. Several excitation lines are resolved due to the 
presence of resonant microwave excitation. We fit the 
five strongest lines with Lorentzian profiles. The ex- 
tracted peak positions are plotted in Fig. |3jb) along 
with a global fit to the model based on Eq. ([T]), 
yielding Stark coefficients Ad\\F\\/Va = 0.11(1) GHz/V 
and d^F^/Va = 0.26(2) GHz/V. These coefficients 
are about four times smaller than those realized un- 
der strong 532 nm excitation, consistent with recent 
observations of enhanced electrical tuning due to pho- 
toionization [15', T6j . 

We note that the average linewidth for single scans 
[38J was Fss = 0.14(3) GHz for NVl and F^s = 
0.48(8) GHz for NV2. In both cases, Fgg is much broader 
than the natural linewidth, Fnat ^ 13 MHz, and is inde- 
pendent of scan rate up to ~ 20 GHz/s. This property 
requires further investigation, as Lgs ^ Fnat has been 
observed elsewhere [14, 16, 38j. 

A likely cause for NV spectral diffusion is charge dy- 
namics due to photoionization of nearby defects. To in- 
vestigate, we use PLE spectroscopy in a different device 
on the |0) \Ex) transition of a single NV center in nat- 
ural, type Ila (Ural) diamond. This sample was chosen 
due to the much narrower linewidth, Fgs = 60(7) MHz, 
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FIG. 3. (a) PLE spectra versus applied voltage for NVl. 
A 532 nm repump pulse (2.5 /xW, 1 s duration) was applied 
every 60 s. (b) Fitted transition frequencies from (a) and 
global fit based on Eq. ([T]), with transitions labeled, (c) PLE 
spectra for a single NV in the Ural sample (|0) \Ex)) at 
Va = —4 V. Green arrows indicate when the repump was 
applied, (d) PLE spectra for the same NV with voltage 
feedback applied, (e) Fitted peak positions for scans in (c) 
[upper panel] and (d) [lower panel]. The recovery under 
feedback after two separate repump pulses are inset, (f) 
Histogram of the peak positions in (e). (g) 100 s segment of 
PLE spectra on NV2 obtained by rapidly scanning a voltage, 
Vac, with feedback applied to Vdc- 



even after ~ 15 MHz of power broadening. Figure [sj^c) 
shows typical PLE spectra for 200 ms scans with re- 
pump pulse (~ 10 /iW, 20 ms duration) applied only 
after the NV center had photoionized. The transition 
frequency drifts over a range significantly larger than 
Tss during the 280 s data set. 



Our solution to the spectral-drift problem is to ac- 
tively adjust Va to compensate for the changing local 
field. We start with Va = —4 V, and, during the back- 
scan of subsequent scans (final 10% of each cycle), we 
employ software-controlled feedback with the following 
algorithm. We first determine the position and inten- 
sity of the peak fluorescence. If the intensity falls be- 
low a threshold, we apply a repump pulse and do not 
change Va- Otherwise, we change Va based on optimized 
proportionality and integration inputs [Supplementary 
Information] . 

Figure |3jd) shows PLE spectra under similar condi- 
tions as Fig. [sjc) but now with feedback applied. While 
Tss remains unchanged, the center-frequency drift is 
substantially reduced. We fit the spectra in Figs, [sj^c) 
and (d) with Lorentzian profiles and plot the extracted 
peak positions in Fig. |3je). In the case of no feedback, 
two mechanisms of spectral drift are identified: large in- 
stantaneous jumps following application of the repump 
and slower drift in between repump pulses. Under feed- 
back, spectral jumps still accompany repump pulses, 
but these are quickly compensated for. Two spectral 
jumps with the slowest recovery, 5-10 scans, are shown 
as insets. A figure of merit for the total drift is obtained 
by plotting the histogram of fitted peak positions from 
all PLE scans and determining the resulting standard 
deviation, a [Fig. |3jf)]. Without feedback, we find a 
non-uniform profile with a ~ 65 MHz. Under feedback, 
a = 16 MHz, which is smaller than Tss and comparable 

to Fnat- 

This feedback technique can be applied at signifi- 
cantly higher bandwidth (here, up to 20 Hz scan repe- 
tition rate) without compromising stability. Through- 
out, we find that feedback reduces <j to a fraction of 
Tss- Similar results were obtained for NVl and NV2 
[Supplemental Information] as well as for stabilizing the 
I ± 1) ^ ^2 transition. 

It is often advantageous to perform experiments with 
the excitation laser frequency fixed to an external ref- 
erence. In this case, voltage feedback can still be em- 
ployed by sweeping the ZPL transition frequency using 
an AC voltage. Vac, and providing stabilizing feedback 
to the DC component, Vdc- With this technique, feed- 
back can be applied continuously without substantially 
degrading photon indistinguishability, provided that the 
modulation depth and laser linewidth are sufficiently 
small. 

Figure [3]^ g) shows results of locking the NV2 |0) 
\Ex) transition frequency using only applied voltages. 
We perform PLE spectroscopy as before except, instead 
of scanning the laser frequency, we ramp the voltage. 
Vac (applied to V2 and Kef) 7 with amplitude 3 Vpp and 
period 0.1 s. Meanwhile, Vdc is fed back to Vi, initially 
starting at —8 V, but varying by ~ ±4 V throughout 
the 600 s measurement. After background subtraction, 
we collect on average 144 cts/s. This compares favor- 
ably to the 34 cts/s collected without feedback (with 
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repump applied every scan). The overall count rate can 
be further increased with improved collection efficiency 
[47l H8\ and resonant Purcell enhancement [28, 49j. 

In summary, we have used the Stark effect to elec- 
trically tune and stabilize the structure of the NV cen- 
ter's excited state. Applied simultaneously to a pair of 
NV centers, these techniques pave the way for increased 
two-photon interference visibilities [TTl [131 US] her- 
alded entanglement success probabilities [8]. 
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Supplementary Information: Dynamic stabilization of the optical resonances of single 

nitrogen-vacancy centers in diamond 



ELECTROSTATIC MODELING 



We modeled the electric field produced by the elec- 
trode structure in Fig. 1(c) of the main text using the 
COMSOL MULTIPHYSICS® electrostatics package. 
The bottom layer of each electrode is composed of 10 
nm of Ti, so for simplicity we model the electrodes as 
being 100-nm thick, composed entirely of Ti. We use a 
relative permittivity for the diamond substrate = 5.1. 
The positions of NVl and NV2 were determined by flu- 
orescence micrographs, as in Fig. 1(d) of the main text. 

Based on our simulations, application of 10 V to 
one electrode (with the other two electrodes grounded) 
corresponds to an electric field amplitude at the lo- 
cation of NVl of 0.9,0.6, and 1.1 MV/m, for Vi,V2, 
and Vref? respectively. In Fig. [ST]we plot the in-plane 
electric field components on the diamond surface for 

{T^l,V^2,Kef} = {10,0,10} V. 

As a note of caution: this model assumes a perfect 
dielectric response. In reality the local field is subject 
to significant deviations due to charge variations intro- 
duced by the electrodes or from photo- ionization [SlJ. 
These deviations do not substantially affect the perfor- 
mance of dynamic ZPL stabilization, since the propor- 
tional gain can be adjusted to compensate (see below), 
but they do play an important role in static tuning. 
One example is the photo-induced effect that is respon- 
sible for the ~ 4 times greater Stark tuning coefficients 
for strong (2 mW) cw green excitation, compared with 
weak (~ 60 nW) red excitation, and is discussed in de- 
tail in Ref. [SlJ. 

We also observe hysteretic charging effects which al- 
ter the static tuning even with very weak (< 60 nW) 
red excitation and without applying any green repump. 
This hysteresis results in a variation of observed tuning 
coefficients which depends on scan direction, rate, and 
history. The variation exists even between successive 
voltage scans as short as 10 s. Overall, the observed 
tuning coefficients vary by up to a factor of 5, making 
a measurement of the intrinsic dipole moment of the 
NV center a difficult task. Using the Stark tuning co- 
efficients for NVl, determined from the 2 hr forward 
voltage scan in Fig. 3(a) of the main text, the dipole 
moment is {/S.d\\^d±_} ~ {4,5} GHz/MV/m, which is 
probably correct to within a factor of 2-3. For this cal- 
culation the magnitude of the applied electric field and 
its angle with respect to the NV axis were estimated 



-2 2 MV/m 




/\ F 

/ NV2 ^ y' 
NVl ^ 




lOjxm Yf^,F,,FJ={10,0,10}V 

FIG. SI. In-plane electric- field components on the diamond 
surface, F^' ^y> ^ for Va — ^0 Y applied simultaneously to Vi 
and Kef, with V2 = 0. The locations of NVl (dashed circle) 
and NV2 (asterisk) are shown, as in Fig. 1 of the main text. 
The out-of-plane component F^' ~ 0. 



REPRODUCIBILITY OF EMISSION 
SPECTROSCOPY 



The reproducibility of Stark emission spectroscopy 
results was tested by varying the voltage on various 
combinations of electrodes. Figure [S2ja) shows emis- 
sion spectra while varying the voltage applied to Vi, 
with V2 = Kef = 0. Despite the continuous voltage 
tuning, small kinks in the emission lines are observed, 
particularly following dark frames. During the dark 
frames of this particular data set, we performed a peak- 
finding procedure that locates the optical focal position 
which produces maximum NV emission. The resulting 
shifts in optical focus may produce small changes in the 
local field due to photo-induced charge redistribution 
|S H IS2 | . In Fig. 2 of the main text we used a different 
procedure to maintain optical focus based on continu- 
ous feedback using a weak white-light reflection image. 
This may explain the absence of sharp kinks in the Stark 
emission spectra in Fig. 2(a). 



The inset of Fig. |S2[a) shows the low-field spectra, 
with three emission lines clearly visible. As in Fig. 2 of 
the main text, we do not observe IA2) ^ | ±1) emission, 
indicating ground-state spin polarization, ^gs 



based on the electrostatic modeling in Fig. [ST] and the 
dependence of NVl fluorescence intensity on excitation 
light polarization. 



> 

Following the procedure outlined in the main text, we 
fit these spectra [Fig. [S2|b)] and found Stark coeffi- 
cients of Ad||F||/K = 0.26(2) GHz/V and d^F^/Va = 
0.81(4) GHz/V. In Fig. [S2|c), the total emission in- 
tensity as a function of 5± is plotted. These data are 
qualitatively similar to that found in Fig. 2(d) of the 
main text. In Fig. |S2[d), the emission intensity of the 
\(oy) I i 1) line, relative to the total emission from 
\(oy), is plotted versus S_\_. Based on the fit, we find 
Or = 12(5)°. 
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FIG. S2. (a) Stark emission spectra as a function of applied voltage on Vi, with V2 = ^ref = 0. Low- field data are displayed 
on the right, (b) Low- field peak positions determined from Lorentzian fits (symbols) and global fit. Lorentzian fit uncertainty 
is smaller than the plotted symbols, (c) Total ZPL emission versus 6±. (d) Relative intensity of the \S'y) ^ | zb 1) emission 
line along with fit. The methods for fitting and processing spectra are described in the main text. 



FEEDBACK OPTIMIZATION 

The basic feedback protocol is described in the text, 
and here we outline experimental details. All PLE scans 
used a ramp waveform with 90% duty cycle. In other 
words, we scan either the laser frequency or Vac in one 
direction for 90% of the total scan cycle and the final 
10% is devoted to scanning back in the other direction 
(the "back-scan"). We divide our collected FSB counts 
into bins of variable size, typically forming 10-50 bins in 
total. During the backscan of each cycle (denoted with 
index z), we search for the bin location, 6^, with the 
maximum counts, Q. We set a threshold, T, typically 
corresponding to a count rate of 1 kcts/s, much larger 
than the background signal off resonance. If Q < T, 
we apply a green repump pulse for the remainder of the 
backscan and do not change Vdc- 

If Q > T, we do not apply a repump pulse. Instead, 
we change Vdc by an amount SVi using the following 
formula: 



1 ^"^ 

6V, = Gx{B--Y,b,_,). 



(SI) 



Here G is a gain factor, N is an integration factor corre- 
sponding to the number of cycles used to determine SVi, 



and B is the desired peak position. In our experiments, 
we typically set B to be the bin at the center of each 
scan. For all of the laser-frequency scans in the main 
text [Fig. 3(d)-(f)], we used N = 1 (no integration). 
For the 0.1-s voltage scans we found that feedback was 
most efficient using TV = 2-4. The 100-second portion 
shown in Fig. 3(g) used N = 2. We separately optimize 
G based on the NV center's Stark tuning coefficients as 
well as the method and rate of scanning. Figure [SSj^a) 
depicts a timing diagram of the feedback routine. 

We performed the feedback routine discussed above 
on the |0) \E^) transition of NV2 in the NV-doped 



surface layer sample. Figure S3^b) shows typical PLE 
spectra when scanning the excitation laser frequency 
through resonance and applying a repump after each 
scan. As mentioned in the main text, the average 
linewidth for single scans, computed using the technique 
in Ref. [S3J, was F^s = 0.48(8) GHz for NV2. Nonethe- 
less, the ZPL center frequency in Fig. [S3jc) drifts over 
a range significantly larger than Fgs during the 300 s 
data set. 

^c) shows PLE spectra under similar con- 



Figure S3 



ditions as Fig. [S3[c) but now with feedback applied. 
While Egg remains unchanged, the center-frequency 
drift is substantially reduced. A figure of merit for 
the spectral drift of the transition is obtained by sum- 
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FIG. S3, (a) Timing diagram of the feedback protocol de- 
scribed in the text, (b) PLE spectra for NV2 (|0) \Ea:)) 
with Va = —4 V. The scan rate was 1 Hz, and a repump 
pulse (10 /iW, 0.1 s) was applied during every back-scan, (c) 
PLE spectra for NV2 with feedback applied. The enhanced 
spectral stability allowed us to reduce the laser-frequency 
scan range by a factor of four, (d) Sum over the scans in 
(b) and (c) with, respectively, feedback off (blue) and on 
(red) along with Gaussian fits (solid black lines) using the 
mean background of 150 cts/s (black dashed line), (e) Com- 
parison of Finhom as a function of scan repetition rate. The 
spectral-diffusion-broadened linewidth (blue shaded region) 
is the range observed over data sets taken without feedback 
with 1, 2, 5, 10 and 20 Hz scan rates. 



ming over many PLE scans and determining the re- 
sulting inhomogenous linewidth, Finhom- This figure-of- 
merit is somewhat different from the histogram tech- 
nique described in the main text and was employed 
due to the large single-scan linewidth for this NV cen- 
ter. Figure [SSjd) compares the sum of spectra with 
feedback off [Fig. [SSj^b)] and on [S3[c)] along with 
Gaussian fits. Without feedback, we find Finhom = 
1.8(2) GHz ^ 3.8Fss, and with feedback this decreases 
to 0.54(4) GHz ^ l.lFss. 

This feedback technique can be applied at signifi- 
cantly higher bandwidth without compromising stabil- 
ity. Figure [SSj^e) plots Finhom as a function of scan 
repetition rate up to 20 Hz. Throughout this range, we 
find 

Tinhom ^ l-3Fss. Futurc improvements could in- 
volve ultra-fast correlation measurements to determine 
the nature of the broad single-scan linewidth [S4] . 
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